On intraband superconducting instability smeared by interband interaction 
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We study several manifestations of former intrinsic phase transition in the band with weaker 
superconductivity by analysing thermodynamics of a two-gap system. We show that weak interband 
interaction acts as external field on that critical point. At that, the corresponding singularities of 
the thermodynamic functions become smeared and all located at different temperatures. 
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I. INTRODUCTION 

After an extension of the Bardeen-Cooper-SchriefFer 
theory to the case of overlapping bands^ multi- 
component superconductivity was confirmed for many 
groups of materials: metals^, doped SrTiOg^, cuprates*, 
MgB2&"~ and pnictides* among others. 

In several two-band systems both gaps behave with 
temperature in a very similar manner throughout the 
whole superconducting phase^'». At the same time, a 
two-band model has distinct intraband critical points. 
Therefore, there exists a temperature region where 
stronger band (i.e. band with stronger superconductiv- 
ity) is intrinsically superconducting (active), but weaker 
one is not yet (passive). Arbitrarily weak interband cou- 
pling fixes single phase transition for joint superconduct- 
ing condensate by penetration of the Cooper pairs from 
an active band into the passive one inducing there super- 
conductivity (interband proximity effect). As a result, 
for tiny interband pairing the gap in the weaker band in- 
flects somewhere below critical point of the joint conden- 
sate and then slowly reduces until the phase transition 
temperature is reached. 

There exists a number of experimental verifications 
for the behavior of that type. Almost vanishing inter- 
band interaction was observed in Nb doped SrTiOg-iiS 
and in NbSe?^^. The conductivity and superfluid density 
data points to the nearly decoupled bands in VsS i^^i^^ 
and FeSei_ar^. The studies of nickel borocarbidesi^ and 
sesquicarbidesi^ are in good agreement with prediction of 
weakly interacting two-band scenario. The dependence of 
resistive superconducting transition on the current direc- 
tion in UNi2Al3 is explained in terms of weakly coupled 
order parameteraii. It has been revealed that the gap 
in the inner planes of multilayered system Cul234 fully 
opens below 117 K, but that in the outer planes only be- 
low 60 K"'^^. An impact of the coupling between order 
parameters on the superconductivity in layered meso- 
scopic rings was illustrated by proximitization between 
two participating layers-^. The change of interband in- 
teraction with doping was demonstrated in Mgi_2;Al2.B2 
leading to the appearance of a "tail" at the Tr-gap^". 
Similar feature was found also in iron-based oxypnictide 

LaFeAsOi^^F^r^^ 

In spite of single superconducting instability of a two- 



band system the memory about intrinsic criticality of 
the weaker band can be still visible in the heat capacity 
curve^^i^, coherence lengthsSi, and relaxation times— 
of order parameters fluctuations as additional maximum 
below phase transition temperature. An impact of tiny 
interaband pairing on the separation between vorticea^^, 
vortex size^, and type- 1.5 superconductivity^^ is also 
noteworthy. 

The evolution of the gap healing length, analysed re- 
cently for two-band scenario^^, reveals that interband 
coupling plays the role of the source field governing in- 
trinsic criticality of the weaker band similarly to an exter- 
nal magnetic field for ferromagnetic materials or an ex- 
ternal electric field for ferroelectric systems. More direct 
characteristics of superconductivity should be examined 
to confirm that important conclusion. 

In the present contribution we perform a detailed anal- 
ysis of the smearing of intraband superconducting insta- 
bility seen in the temperature behaviour of thermody- 
namic functions. In this respect the relevant two-band 
materials are discussed. 



II. THERMODYNAMICS OF A TWO-GAP 
SYSTEM 



We start with the Hamiltonian of a homogeneous 
two-band superconductor with intra- and interband pair 
transfer interactions, 

H = ea(k)a+^^aQks 
7 W^aa'aakfaa-k^a^'-k'-Laa'k't, (1) 



kk' 



where = — /i is the electron energy in the band a = 
1, 2 relative to the chemical potential ijl;V\s the volume 
of superconductor and Waa' are the matrix elements of 
intraband attraction {Waa > 0) or interband (a ^ a') 
interaction. It is supposed that the chemical potential 
is located in the region of the bands overlapping. We 
assume that (effective) electron-electron interactions are 
nonzero only in the layer ^,±hwD and that the interaction 
constants are independent on electron wave vector in this 
layer. We also take W12 = W2i- 
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To calculate the partition function Z — Spexp i^^-^^p^ 

we use the Hubbard-Stratonovich transformation^ by 
means of which one linearizes and diagonalizes an ex- 
ponent in Z by introducing complex integration vari- 
ables. For sufficiently small interband interaction = 
W11W22 — > the static path approximation reads 



Z = 



V 



exp 



Vf 



TTknTW J J " \ kBT 

— 00 ^ ^ 

^^fg„(k)-2fcsTlnf2ch^"^^) 



dS^dd'IdS'^dS'l (2) 



2kBT 



(3) 



Here the integration variables S'^ and S'^ are treated as 
real and imaginary parts of the non-equilibrium complex 
order parameters Sa, and 



?a(k) 



2V 



\Sc.\ 



(4) 



Next we find the equilibrium free energy density / = 
— k^T^^ for a macroscopic superconductor {V -> 00). 
First, we perform an integration in Z over the phases 
of non-equilibrium order parameters and then go to infi- 
nite volume. In this process one obtains the integration 
over wave vector k instead of summation in /, which 
we replace with the integration over energy. Free energy 
density / becomes 



Ak-eTpa / In 



ch 



2kuT 



ch 



den 



-A^ TTTT^ AaAa-a 



(5) 



A^, where Aq is the modulus of 



Here Ea{x) — \/ x 
equilibrium order parameter, i.e. the value of I^qI which 
minimizes non-equilibrium free energy density /. The 
quantity /„ corresponds to the free energy density in the 
absence of superconductivity, and pa is the density of 
electron states taken to be constant in the narrow inte- 
gration layer around the Fermi level. 

The minimization of non-equilibrium free energy leads 
to the equations for equilibrium order parameters 



3— a, 3— a 



Aq.-/9qAc 



th 



-Ea(g) dg 

2kBT Eo,{i) 



3-q| 
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If interband interaction is absent, the latter system splits 
into two independent equations which describe intrinsic 
superconductivity in the bands. The corresponding crit- 



ical temperatures equal Tc, 



1.13 



How- 



ever, if interband coupling is present the critical temper- 
atures Tea transform into 



1.13— — e 



(7) 

where Tc- > Tc+. With \Wi2\ increase the tempera- 
ture Tc- increases and Tc+ decreases approaching zero 
as 0. For W12 one has Tc- Td and 

Tc+ — ^ Tc2, if Tci > Tc2- Note that for a superconductor 
with interacting bands there is only one phase transition 
temperature Tc = Tc-. However, the point Tc-\- is also 
meaningful, because below Tc-\. the metastable supercon- 
ducting states or saddle-points of non-equilibrium free 
energy appear—. These peculiarities can be reflected in 
the behaviour of superconducting fluctuations. 

Thermodynamics of the system related to supercon- 
ductivity is entirely described by excess free energy den- 
sity A/ = f — fn and its temperature derivatives, e.g. 
excess specific entropy As = —A/' and specific heat ca- 
pacity Ac = —TAf". For the degenerate electron gas 
one has Sn = c„ = JsT, where 73 = |7r^A:|(pi -I- P2) is 
the Sommerfeld constant. Therefore, the extrema and 
inflection points of As coincide in the temperature scale 
with zeros and extrema of ^ , correspondingly. It is also 
worth to mention that the quantities A/, As and Ac 
may be represented by means of additive contributions 
from the relevant bands. 



III. RESULTS AND DISCUSSIONS 
A. Numeric study 

Below we examine thermodynamic properties of a two- 
gap superconductor with extremely weak interband cou- 
pling. In such a system the joint superconducting state 
is close to the splitting into two independent conden- 
sates. For illustration we consider the following set 
of intraband parameters: pi, 2 = (l,0.9)(eV • cell)~^, 
W^ii,22 = 0.3 eV • cell and cell — 0.1 nm'^. In this case 
Tc2 = OMTci. 

Fig. [T] shows the changes of the gaps, excess entropy 
and heat capacity with intra- and interband interactions 
calculated numerically on the basis of Eqs. ([S])-®. By 
turning interband coupling on, the intrinsic criticality in 
the weaker band becomes smeared by interband proxim- 
ity effect and both gaps go to zero at very same point 
Tc = Tc-, i.e. the vanishing of the smaller gap A2 at Tc2 
changes into inflection near Tc+ ~ Tc2- That inflection 
indicates the crossover from active to passive regime of 
the weaker band and it disappears as interband interac- 
tion exceeds some value. 

To interpret the temperature behaviour of entropy and 
heat capacity it is convenient to visualize band contri- 
butions, see panels 2b, c, and 3b, c in Fig. [TJ For tiny 
interband coupling these contributions start to develop 
considerably below different temperatures ~ Tci^2 
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FIG. 1. The plots of the gaps (a), excess specific entropy (b) and normalized excess specific heat capacity (c) vs reduced 
temperature for 1^12 = 0.0001 eV • cell (1), W12 = 0.001 eV ■ ceU (2) and W12 = 0.01 eV ■ cell (3). In the first row of figures 
the dashed lines denote the corresponding quantities for slightly weaker coupling in the second band W22 = 0.25 eV • cell 
(Tc2 = O.SSTci). In Fig. (la) the larger gaps flow together for the constants 14^22 considered. In the second and third rows of 
figures the thin lines denote band contributions to the excess entropy and specific heat capacity. 



resulting in the substantial increment for entropy and 
the non-monotonicity for heat capacity in the vicinity 
of Tc+. These features demonstrate again the critical- 
ity of the weaker band smeared by interband interac- 
tion. The latter was observed also experimentally^. For 
sufficiently strong interband pairings both bands have 
comparable contributions to the joint superconductiv- 
ity up to Tc- Note that the presence of several compo- 
nents of superconductivity makes the standard Bardeen- 
Cooper-Schrieffer theory inappropriate for the systems 
considered^i^. Well-known universal ratios become also 
violated22i2i. 

Interestingly, the certain configuration of intraband in- 
teraction channels together with sufficiently weak inter- 
band coupling result in the local maximum of entropy in 
superconducting state, see the first row of panels in Fig. 
[TJ In this case ^ ^ crosses zero below Tc three 
times. Similar effect was revealed in Cul234^. 

Whereas both sharp non-monotonic behaviour of heat 
capacity and the kink of the smaller gap are the foot- 
prints of intrinsic phase transition in the weaker band, 
should one identify the maximum of ^ below Tc with 

Tj^f , the temperature where A2 inflects? Numeric cal- 
culations demonstrate (see Fig. [2^) that these points 
match only approximately and not for extremely weak 
interband interactions. The minimum of — below Tr is 

more farther off T,^/. We see that neither the inflections 
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FIG. 2. (a) Nearest to Tc+ inflection points for A/' (two 
thin dashed curves for maximum and minimum of ■^), A/" 

(thin solid curve), A/™ (thick solid curve) and for A2 (thick 
dashed curve) vs interband constant, (b) Inflection points for 
A2 (thick dashed curve) and for A/"' (thick solid curve) com- 
pared to the approximations given by the analytic formulas 
(fTTjl and (O (thin solid curves). 



for A/'(T) nor for Af"{T) reproduce the position of the 
kink of A2(r), however, it is the function Af^{T) which 
inflects at T^^f for extremely weak interband pairings. In 
the vicinity of Tc+ the contribution from stronger band 
to A/^(r) is practically temperature independent and 
remarkably smaller than that from weaker band. The 
latter function has a maximum at the temperatures con- 
sidered which flattens as interband interaction increases. 
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Material 


/ 12 




Ac 

c„ _^A2 


Ac 

c„ _„A2 
-'min -'inf 


„A/" „A2 
-'int "-'int 


yA/III_j,A2 
int inf 


V W11W22 


Tc 






Tc 


LiFeAs 


0.096=^^ 


0.310 








-0.175 


VaSi 


0.014iS 


0.490 


0.009 


0.102 


0.050 


-0.021 


FeSei-, 


0.003i^ 


0.376 


-0.019 


0.062 


0.017 


-0.004 



TABLE I. The relative position of the nearest to Tc+ inflection points for A/' (both majdmum and minimum of ^ are 
considered), A/" and for A/"" in case of LiFeAs, VaSi and FeSei-^,. Material parameters used in numeric calculation -were 
taken from the references shown. 



Next we compare the positions of the relevant pecu- 
liarities for two-band superconducting materials under 
interest, see Table ID In LiFeAs heat capacity behaves 
monotonically below and T-^^ can be found only ap- 
proximately, i.e. in this system W12 is sufficiently large. 
VaSi can be considered as a compound with moderately 
weak interband coupling due to the closeness of T-^f and 
the position of maximum of ^ . This type of behaviour 
is described by r.h.s. of Fig. FeSei_^ is a system 
where the kink of the smaller gap appears close to the 
inflection point of /S.f^ (middle part of Fig. [2^). Our 
estimations indicate also that entropy of FeSei-a: should 
have a local maximum in superconducting phase approx- 
imately at 0.41Tc = 3.4 K. This feature was discussed 
above in connection with dashed lines in the first row of 
panels in Fig. [1] 



B. Analytic treatment 

Below we confirm numeric result about the discrep- 
ancy between the positions of the kink of the smaller gap 
and heat capacity maximum in the case which can be 
treated analytically. By starting with non-equilibrium 
free energy / expanded in powers of 5a, we obtain for 
macroscopic superconductor 

M = Y. («"^' + y - 7A.A3_.) , (8) 



where 



and ba 



W- 



3 — a, 3 — a 



1^2 



- Pa hi 
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(9) 



O.llpa 
(feBT)^' 
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\Wl2\ 



This approximation works 
perfectly in the vicinity of critical point. However, as we 
do not expand the parameters Oq and ha in powers of 
Tc — T, we have at least qualitatively correct picture also 
further off the critical temperature^. In particular, one 
ascertains numerically that corresponding gap equations 



(10) 



still describe an inflection in the behaviour of A2(T), 
while the derivatives of expressions (|8]) and ([5]) are quite 
close for relevant temperatures. Note that alternative 
approaches based on the expansion in powers of — T 
cannot reproduce the kink of the smaller gap2^. 



Next we solve analytically the system (|10p in the vicin- 
ity of Tc+ for extremely small interband couplings, see 
Appendix. The solution found for A2(T) is character- 
ized by the following inflection point 



9y 



35 



{1} 



2A 



{0}2 



(11) 



where 



= 9.4(fcBTe2) 



(12) 



W12 



p^w^^w^i I ■ T^*^ temperature T.^l 



and i = In 2ii , w = { - 

contains the lowest order correction to Tc_|_ proportional 

4 

to Wi2- It proves that the smaller gap as a function of 
temperature inflects in the vicinity of Tc+ if interband 
interaction is extremely weak. 

By using solutions (T) we analyse the temperature 
dependence of free energy ([5]). For the inflection point of 
Af^{T) we get 



^inf 



Tr- 



1 



3B 



{1} 



(13) 



Analytic expressions (jlip and (jl3p approximate values 
found numerically strikingly well, see Fig. [2}d. Our calcu- 
lation confirms also that functions A/™(T) and Af^{T) 
vanish remarkably farther off Tc+ as opposed to A/^(T). 
In other words, the infiection points of entropy As or 
heat capacity — do not appear sufficiently close to 
for tiny interband interactions 

The temperatures Tj 
([T^ coincide, if 



c+ 



A2 



and given in Eqs. (fTTj)- 



3y 



{^}3 



3t < 1, 



(14) 



and this condition is fulfilled when Tci^2 are close. Such 
a two-gap system demonstrates analytically the smaller 
gap and the third derivative of free energy inflecting both 
at the very same temperature as interband coupling be- 
comes vanishing. On the other hand, one obtains 



T.; 



inf 



t: 



■■inf 



(15) 
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i.e. an impact of the discrepancy between temperatures 
7ci,2 on the closeness of inflection points under interest 
may be compensated by the weakness of interband in- 
teraction. Our numerics presented in Fig. [5] verify that 
conclusion. 

Last but not least, in the lowest order in w we have 
c2)^=2/y^i.e. 



A2(r, 
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(16) 



Appendix: Solution in the vicinity of Tc+ 



Therefore, the "susceptibility" related to the "field" W12 
reads as 



dWi2 



(17) 



Moreover, the lowest order term in A/'(Tc2) stemming 
from the weaker band is 



2 P2 2 



c2 



(18) 



The corresponding term in A/"(Tc2) related to the same 
band equals 



2 P2 .{0} 



W^2- 



(19) 



In the Landau mean-field theory of criticality the order 
parameter changes as /i^ {S = 3), susceptibility as h'^~^, 
entropy as /i^ (e = |) and heat capacity as /i~"'= {ac = 0) 
under applied external field h. Thus, interband coupling 
plays a role of external field governing the intrinsic criti- 
cality of weaker band. This consequence is in agreement 
with the critical exponent found for the corresponding 
healing length^^ . Note that intrinsic superconducting in- 
stability of the stronger band is affected in a qualitatively 

different way, for instance, Ai(Tci) ^ 



1 



^ - W^12 



for vanishing interband interaction. That pairing rather 
shifts intrinsic phase transition from Td to Tc, but not 
smears it as in the case of weaker band. 



IV. CONCLUSIONS 

We have studied the interrelations between the ther- 
modynamic features of two-gap superconductors (e.g. 
FeSei^x, VaSi) caused by former intrinsic phase tran- 
sition in the weaker band. In particular, there appear 
atypical maxima in the temperature dependencies of en- 
tropy and heat capacity, and the gap kink in the weaker 
band. These peculiarities have slightly different locations 
on the temperature scale, and can move away from each 
other with decrease of interband coupling. The latter in- 
teraction is shown to play a role of external field which 
smears intrinsic superconducting instability of the weaker 
band. 



By introducing the notations x 



Al A 



-fa and 



— — one gets from system (fTO|) the following equa- 



y 

tion 



;((x-A)3 



0102) 



(A.l) 



As the condition 7^ — 0102 = determines the points 
Tc±, we obtain ~ A+ + ?/-|_ (index + implies T = 



V°° A 



{!+»} 



and 



Next we denote A^ 

X]^o 2/+^"''"'^' ^'^'^ summation is made over the contri- 
butions of different powers of w (upper indexes in curly 
brackets). We also introduce the quantities A ^ ^ 



for which Aa 



and B = ^ 

EDO + 

given in the main text, and 



{«} 



and Bs 



The first terms of these expansions are 



.{1} 



Pi 



w. 



(A.2) 



By differentiating Eq. (jA.l[) at Tc+, we find the con- 
secutive temperature derivatives of x. For instance, in 
the lowest order in w 



1 -24">yy^'+62/;^ 



{3}3 _^{0}^{l} 



c+ 



2 4{0}2 



.{1} 



rp2 
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,(A.3) 



(A.4) 



etc. From these expressions one gets the value of A2 and 
its derivatives at . The analytic dependence A2 (T) is 

given then by Taylor series A2+ -|- A'2+{T — Tc+) + 

Usually w power of the lowest order contributions to 
the consecutive derivatives at Tc+ decreases by ^ as or- 
der of the derivative increases, for instance, ^ ^ wi or 

A™ 1 

-rj^ ^ w -i . However, sometimes corresponding w power 
decreases by | due to extra vanishing of the terms. One 



notable example here is 



A2+ 

A, 



w 3 , resulting in the solu- 



tion of the equation A2{T^^f) = in the form 



t; 



inf 



Tc+- 



A" 

^2+ 

AH 
^2-1- 



Tc+\ 1 



3B 



{1} 



2A 



{0}2 



2Af.?^ 



(A.5) 
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This value contains the lowest order correction to Tc+ 

2 i 

proportional to W3 ^ Wy2- Note that zero (or negative) 
power of that correction disagrees with an assumption of 
the closeness to Tc+ in case of small Wi2. 

By using free energy ([8]) rewritten in the form 

A/ = 62x(-x+^ + |^(a;-Af), (A.6) 

we find that the situation with the consecutive temper- 
ature derivatives of A/ at Tc+ is similar to that de- 
scribed above for the derivatives of the smaller gap. Here 
the suitable ratio of the derivatives occurs only for their 

A/V 2 

higher orders, . ,vi ^ w a , resulting in the inflection point 



of A/™(r) in the form 

i.e. the lowest order derivative of Af{T), which goes to 
zero sufficiently close to Tc+, is A/^(T). 

Finally, for weak interband coupling one has ?^ 1 + 

oiw). Thus, the values A2(rc2) = EZo ^^^+ ^^"^^^ 
or A/M(re2) = A/|"+™)(«±i:: (here n and 

m are the orders of the temperature derivative) con- 
tain the sum of w dependent contributions which be- 
come smaller and smaller as n increases. For tiny W12 
these values are defined predominantly by A2(T'c+) or 
A/('")(Tc+), correspondingly. 
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